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Abstract. It has recently been suggested by Luminet et al. (2003) that the WMAP data are better matched by a 
geometry in which the topology is that of a Poincare dodecahedral model and the curvature is "slightly" spherical, 
rather than by an (effectively) infinite flat model. A general back-to-back matched circles analysis by Cornish et al. 
£C) ' (2004) for angular radii in the range 25 — 90°, using a correlation statistic for signal detection, failed to support 

this. In this paper, a matched circles analysis specifically designed to detect dodecahedral patterns of matched 
circles is performed over angular radii in the range 1 — 40° on the one-year WMAP data. Signal detection is 
attempted via a correlation statistic and an rms difference statistic. Extreme value distributions of these statistics 
are calculated for one orientation of the 36° 'screw motion' (Clifford translation) when matching circles, for the 
opposite screw motion, and for a zero (unphysical) rotation. The most correlated circles appear for circle radii 
of a = 11 ± 1°, for the left-handed screw motion, but not for the right-handed one, nor for the zero rotation. 
f*—} ' The favoured six dodecahedral face centres in galactic coordinates are (l U ,b U ) » (252°, +65°), (51°, +51°), 

(144°, +38°), (207°, +10°), (271°, +3°), (332°, +25°) and their opposites. The six pairs of circles independently 
each favour a circle angular radius of 11 + 1°. The temperature fluctuations along the matched circles are plotted 
and are clearly highly correlated. Whether or not these six circle pairs centred on dodecahedral faces match via 
a 36° rotation only due to unexpected statistical properties of the WMAP ILC map, or whether they match due 
to global geometry, it is clear that the WMAP ILC map has some unusual statistical properties which mimic a 
O ' potentially interesting cosmological signal. 
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1. Introduction (based on distributions of gravitationally collapsed objects 

jj] ■ distributed in three-dimensional comoving space). For a 

C3 In the past twelve months, several authors have analysed i ist and discussion of both two-dimensional and three- 

the possibility that the primordial temperature fiuctua- dimensional methods, see Table 2 of Luminet & Roukema 

tions in the cosmic microwave background, as measured ( lg9 g) and the accompanying discussion, 
in the first-year data of WMAP (Wilkinson Microwave 

... -n , n , m i i i nnno j The WMAP data has motivated many two- 

Amsotropy Probe) satellite (Spergel et al. 2003, and ac- J 

\ , , , , dimensional analyses, 

companymg papers), could be better matched by a per- J 

turbed Friedmann-Lemaitre-Robertson- Walker (FLRW) Some of the authors mentioning either multiply con- 
model in which the global geometry is multiply connected nected models consistent with the WMAP data, or indi- 
rather than simply connected. rect hints of multiple connectedness, include Spergel et al. 

The possibility of the Universe being multiply con- ( 2003 ); Tegmark, de Oliveira-Costa, & Hamilton (2003); 

nected is first known to have been suggested by Chian S et aL ( 2003 ); Go P^ Huterer, & Starkman (2003), 

Schwarzschild (1900, 1998). For recent reviews on cosmic " nndm S low values of low 1 multipoles or applying a 

topology, see Lachieze-Rey & Luminet (1995), Luminet multipole vector analysis. 

(1998), Starkman (1998) and Luminet & Roukema (1999). However, other authors (e.g. Feng & Zhang 2003; 

For workshop proceedings on the subject, see Starkman Cline, Crotty, & Lesgourgues 2003; Contaldi et al. 2003; 

(1998) and following articles, and Blanlceil & Roukema Efstathiou 2003a) have suggested various non-topological 

(2000). Detection strategies include both two-dimensional explanations for the low multipole WMAP Ci spectrum, 

methods (based on temperature fluctuations in the sur- such as double inflation or other phenomena from early 

face of last scattering) and three-dimensional methods universe physics, or positive curvature. 
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Yet others (Naselsky, Doroshkevich, & Verkhodanov 
2003, 2004) point out that at least the "internal linear 
combination" (ILC) map of the WMAP data contains 
non-Poissonian signal, due to foreground residues on large 
scales. On smaller scales, Giommi & Colafrancesco (2003) 
find that somewhere between 20-100% of the signal at 
spherical harmonic / values in the range 500 < I < 800 
may be due to blazars. 

Attempts to exclude classes of global geometry mod- 
els using the WMAP data include the calculations of de 
Oliveira-Costa et al. (2004) and of Cornish et al. (2004). 
The latter performed a general back-to-back matched cir- 
cles analysis for angular radii in the range 25 — 90°, using 
a correlation statistic for signal detection. They failed to 
find matched circles for a wide class of models, including 
the torus models, up to a scale of 16.8/i _1 Gpc, which 
Cornish et al. 2004 state as 24 Gpc, since they adopt 
H = TOkms^Mpc- 11 . 

Disagreement also exists on whether or not the low 
quadrupole is really significant in rejecting the infinite flat 
'concordance' model (Efstathiou 2003b, 2004). 

Possibly one of the strongest claims in favour of a pos- 
sible detection is that by Luminet et al. (2003), who point 
out that given standard assumptions on the statistics of 
the fluctuations, the Poincare dodecahedral model implies 
a quadrupole and an octupole very close to those calcu- 
lated from the WMAP data, for the same value of the 
total density parameter, Ot t ~ 1-013 ± 0.02. 

In principle, this is excluded by the Cornish et al. 
(2004) analysis. 

The Poincare dodecahedral model requires positive 
(spherical) curvature. Luminet et al. (2003) favour a to- 
tal density parameter of S^tot ~ 1-013 ± 0.002 based 
on the spherical harmonic statistical analyses of the 
WMAP data, with non-relativistic matter density pa- 
rameter, fi m = 0.28 and cosmological constant 51a = 
f^tot — ^m- This implies that points on the surface of 
last scattering which are multiple topological images of 
single physical points in space-time should correspond 
to matched circles (Cornish, Spergel, & Starkman 1996; 
Cornish, Spergel, & Starkman 1998) which subtend an- 
gular radii at the observer of 'about 35°' (Luminet et al. 
2003). 

Since the angular relations between face centres 
are identical for Euclidean and spherical dodecahedra, 
Euclidean calculations of the relative positions of circles 
are sufficient for testing the Poincare dodecahedral hy- 
pothesis. For example, the Euclidean half-angle of about 
31.7° is valid for angular separations of adjacent face cen- 
tres in the spherical dodecahedron. Since the adjacent cir- 
cles expected by Luminet et al. (2003) have larger angular 
radii (subtended at the observer), this implies that they 
intersect with the face edges. 

However, as is shown in Sect. 13.21 this angular radius 
is extremely sensitive to the value of i^tot- F° r example, 

1 The Hubble constant is parametrised as h = 
//o/100km s" 1 Mpc" 1 . 



keeping fl m = 0.28 fixed, it is sufficient to decrease fitot 
to fitot = 1-009 to bring the angular radius to nearly zero. 

This implies that a Poincare dodecahedral signal may 
have been missed by Cornish et al. (2004) because they did 
not explore the part of parameter space for small matched 
circles. 

In this paper, the missing part of parameter space is 
investigated. A matched circles analysis specifically de- 
signed to detect dodecahedral patterns of matched circles 
is performed over angular radii in the range 1 — 40° on 
Internal Linear Combination map (ILC) of the one-year 
WMAP data. The WMAP data are briefly discussed in 
Sect. [21 While the ILC is unlikely to be ideal for the stud- 
ies of Ci statistics, and some authors (cited above) claim 
correlations with foregrounds, it is hard to see how any 
signal mimicking matched circles oriented in a dodecahe- 
dral pattern could be imposed, either by the construction 
method of the ILC or by foregrounds. 

Signal detection is attempted via a correlation statis- 
tic and an rms difference statistic, similarly to Roukema 
(2000b, 2000a). Extreme value distributions of these 
statistics between a right-handed rotation when match- 
ing circles, a left-handed rotation, and a zero rotation. A 
genuine signal should appear for either the right-handed 
or left-handed rotation, but not both, and should not ap- 
pear for the zero rotation. 

The relation between circle angular radius (a) and lo- 
cal cosmological parameters (f2 m 7^tot) an d the statistics 
used are presented in Sect. [21 

The tentative detection of dodecahedrally distributed 
matched circles with a « 10° and analysis of their statis- 
tical significance are presented in Sect.0J Further discus- 
sions and conclusions are made in Sect.|S] 

For reviews on cosmological topology, see Lachieze-Rey 
& Luminet (1995), Luminet (1998), Starkman (1998) and 
Luminet & Roukema (1999). For workshop proceedings 
on the subject, see Starkman (1998) and the following 
articles, and Blanked & Roukema (2000). For a list and 
discussion of both two-dimensional and three-dimensional 
methods, see Table 2 of Luminet & Roukema (1999) 
and the accompanying discussion. The reader should 
be reminded that while microwave background data is 
still the most popular for topology analyses, consider- 
able work in three-dimensional methods has been carried 
out, including, e.g., Lehoucq, Lachieze-Rey, & Luminet 
(1996); Roukema (1996); Fagundes & Gausmann (1998b); 
Roukema & Edge (1997); Roukema & Blanked (1998); 
Gomero, Reboucas, & Teixeira (2002); Lehoucq, Luminet, 
& Uzan (1999); Fagundes & Gausmann (1998a); Uzan, 
Lehoucq, & Luminet (1999); Fagundes & Gausmann 
(1999); Gomero, Reboucas, & Teixeira (2000, 2001). 

For background on spherical multiply connected 
spaces, apart from the recent analysis by Luminet et al. 
(2003), see Gausmann et al. (2001); Lehoucq et al. (2002); 
Riazuelo et al. (2003) for extremely thorough, in-depth 
mathematical background directly related to the cosmo- 
logical context. For a general background on geometry 
and topology, see e.g. Weeks (2001). 
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Fig. 1. Dependence of the angular size a of matched cir- 
cle radii on the total density parameter f\ot anc ^ the non - 
relativistic matter density parameter S7m for the Poincare do- 
decahedral hypothesis. Clearly only slight changes in £\ t are 
needed for the circle radius to change significantly. The val- 
ues fitot = 1-013, a = 35° expected by Luminet et al. (2003) 
are shown by a star. The word 'approx' is used in the label to 
remind the reader that this is not an exact calculation, since 
Euclidean geometry is assumed. 




a (degrees, approx) 

Fig. 2. An inversion of Fig.Q focussing on the angles found 
to be of interest in the WMAP ILC map. Independently of 
S7m, the slope in the range 5° < a < 15° is dQ^ Q ^/da « 
6xl0~ 5 deg- 1 . 



Comoving coordinates are used when discussing dis- 
tances (i.e. 'proper distances', Weinberg 1972, equivalent 
to 'conformal time' if c = 1). 

2. Observations 

See Spergel et al. (2003) and accompanying papers for an 
introduction to the WMAP data. 
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Fig. 3. Dependence of Qm on a, similarly to Fig.|U The slope 
in the range 5° < a < 15° is dQ m /da « 2.2xl0~ 3 deg -1 , with 
only weak dependence on £\ Q t over this interval. 



While some authors make their own combinations of 
the WMAP maps in individual frequencies, it is hard to 
see how any signal mimicking matched circles oriented in 
a dodecahedral pattern could be imposed, by the method 
according to which the Internal Linear Combination (ILC) 
map is constructed. So this is the data set used here. 

However, in order to obtain conservative results, all 
points within 2° of the galactic plane or 20° of the galactic 
centre are removed from the analysis. 

Moreover, a smoothing of 1° HWHM is applied for 
the main calculation, since if a signal is present, it should 
be present both on super-degree scales and on sub-degree 
scales, but should be easier to detect on the larger scales, 
because the search calculations can be made much faster, 
because the naive Sachs- Wolfe dominates over the Doppler 
effect on larger scales, and because there is likely to be less 
influence from subtle effects which need to be taken into 
account as correction factors, e.g. the special relativistic 
correction for the observer's motion with respect to the mi- 
crowave background, Barrow & Levin (2001); Uzan et al. 



(2002). 



3. Method 

3.1. The Identified Circles Principle 

The identified circles principle was first published by 
Cornish, Spergel & Starkman ( 1996, 1998). This defines 
the set of multiply topologically imaged points for a given 
manifold. 

This set can be generated by considering copies of the 
observer in the covering space placed at distances less than 
the horizon diameter from the observer. The intersection 
of the two surfaces of last scattering (spheres) of the ob- 
server and a copy of the observer is a circle. Since the copy 
of the observer is physically identical to the observer, what 
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appears to be two observers looking at one circle is equiv- 
alent to one observer looking at two circles. 

Hence, each pair of virtual copies of the observer in 
the covering space (separated by less than the diameter of 
the surface of last scattering) implies a pair of circles on 
the surface of last scattering which correspond to identical 
points in space-time. 

3.2. Relating circle angular radius a to tt m , Qtot 

As mentioned above, Euclidean geometry is sufficient 
for calculating the positions of matched circles for the 
Poincare dodecahedral hypothesis of Luminet ct al. 
(2003), even though the model is for positive curvature. 
Errors would occur only if, for example, the angular ra- 
dius of a circle, a, is converted into a comoving distance 
in megaparsecs, or if the ratio of the inverse cosine of the 
ratio of the in-radius, r_, to the distance to the surface of 
last scattering, rgm, is used to infer the angular radius of 
a matched circle, i.e. 

a w cos -1 [ — — ) (1) 
V^SLS/ 

Since the values of f2tot °f interest are very close to 
unity, either of these calculations should only lead to small 
errors. 

In fact, a much bigger uncertainty is due to the sensi- 
tivity of a to the values of fitot ano - ^m- 

The curvature radius can be defined in a strict sense, 
without taking an absolute value of the curvature, as 

This has the advantage of showing that the transition 
from 'slightly' hyperbolic space to flat space to spherical 
space, i.e. from f^ofc < 1 to Sltot = 1 to f2tot > 1 starts 
from imaginary values of the curvature radius, increases 
to larger and larger imaginary curvature radii without 
bounds, goes through a singularity, and decreases from 
unbounded large, real values to smaller and smaller real 
values as fitot increases through positive values. 

It is therefore unsurprising that other physical quan- 
tities related to the curvature radius may show rapid 
changes near the flatness limit. 

To illustrate the relation between a, fitot anc i fW> the 
value of r_, the in-radius of the fundamental domain of 
the Poincare dodecahedron, expressed as 

r_ = OMRc (3) 

is used, along with eq. JQ) and the comoving distance to 
the surface of last scattering 

r SLS = r(n m> fitot.tfo, * = HOO). (4) 

Fig. ^ shows how a few tenths of a percent drop in 
f^tot are sufficient to increase the curvature radius so that 
the in-radius of the fundamental dodecahedron becomes 



as large as the surface of last scattering, reducing the circle 
size from ss 35° to zero. 

So, while Luminet et al. (2003) 's statistical estimates 
of Ci values do suggest that f^tot ~ 1-013 and a w 35°, it 
is clear that smaller f2tot values and circle sizes should be 
investigated. 

In Figs |3 and the dependence of fitot and f2 m on a 
in the angular range found to be of interest in the WMAP 
data is illustrated. Given the fact that the Universe is 
anisotropic (we clearly do not live in a true FLRW uni- 
verse; we live in a perturbed FLRW universe) some vari- 
ation in the size of circles could be expected due to per- 
turbations in the values of f2 m and f^tot- 

The variation shown in Fig. 

dft m /da re 2.2xl0~ 3 deg" 1 , (5) 

is clearly much larger than temperature fluctuations on 
scales of about a radian (Z ~ 2 — 3). 
However, the variation in Otcrt) 

dfWda«6xHT 5 deg" 1 , (6) 

is only slightly larger than that indicated by temperature 
fluctuations. Moreover, the unknown physics of the cos- 
mological constant (or quintessence) makes it hard to a 
priori exclude any fluctuations in J7a, especially with an 
amplitude so close to that of matter fluctuations. 

Hence, a variation in circle size of about one degree 
would be consistent with fluctuations in on scales of 
about a radian of amplitude about 6xl0 -5 . 

Along with restrictions in practical computing power, 
this motivates a search strategy and smoothing on a scale 
of about one degree. 

3.3. Parameter space 

The space of possible orientations of a dodecahedron (cen- 
tred on the observer) and possible circle matches is defined 
by five parameters: three for the orientation of the dodec- 
ahedron, one for the mapping between faces, and one for 
the circle radius. 

The orientation of the dodecahedron depends on three 
continuous free parameters. The definition used here in- 
cludes the galactic longitude and latitude (Z^, fr^) of one 
face centre, and a third parameter allowing a rotation of 
up to 9 = 2tt/5 = 72° around the axis defined by the first 
centre and its opposite. The zero point of the rotation 9 
is arbitrary. 

A given triple (Z , b , 6) can equivalently be writ- 
ten as a quadruple (Zp, ftp, Zp , 6p) representing two non- 
opposite face centres, which uniquely define the same ori- 
entation. 

The mapping of one copy of the fundamental domain 
to the next, or what is most interesting in the present 
case, the matching of one face of the dodecahedron to 
its matching face, has one free binary parameter, 7, the 
handed-ness of the screw motion of the mapping (which is 
a Clifford translation). For a given triple (Z n ,& n ,60, the 
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screw motion can either be left-handed or right-handed. 
The amplitude of the screw motion is it/ 5 — 36°, equal to 
the in-diameter of the fundamental domain in units of the 
curvature radius. 

For a rotation of "fir/5, the only physically valid values 
are 7 = ±1. 

The remaining parameter is that discussed above, the 
angular radius of a circle, a. 

Putting these together, these can be written as five- 
dimensional parameter space: 



{l ll ,b ll ,9 n ,a), 



(7) 



where four parameters are continuous variables and one of 
the parameters, 7, has only two physically possible values, 
±1. 

The binary parameter provides an easy method for 
providing a control test which should demonstrate the 
likely distribution of values of statistical measures of circle 
identity for false matches. 

By setting the Clifford translation as a translation plus 
a zero rotation, i.e. setting 7 = 0, rather than a translation 
plus a ±36° rotation, an extra 4-plane in parameter space, 
in which there is sure to be no genuine signal, is provided. 

So, this extra 4-plane is tested here. 

Strictly speaking, values of (2^,6^) only need to be 
searched over one-twelth of the sphere, since outside of one 
spherical pentagon of the spherical dodecahedron, redun- 
dant points would be searched. For simplicity, the search 
is made for ir > +50°, so that some redundant points 
are expected. 

The full parameter search space is: 

0° < Z 11 < 360° 
50° < b 11 < 90° 
0° < 9 < 72° 
7 G {0,-1, +1} 
1° < a < 40°, 



(8) 



where the zero screw motion is a control experiment which 
can only give false detections. 

3.4. Correlation and difference statistics 

For a given triple (I ,b ,6), a, and the handed-ness 
of the screw motion, there are six pairs of circles which 
should have identical temperature fluctuations if there 
were only the naive Sachs- Wolfe effect, no foregrounds and 
no other statistical or systematic sources of noise. 

Here, the statistic of Cornish et al. (1998) is used as 
the correlator, 



S 



' 5T 



(9) 



where i and j are the hypothetically multiply imaged loca- 
tions in space-time. The mean is calculated using linearly 



interpolated temperature fluctuations in the ILC at inter- 
vals of 0.5 great circle degrees. 

This is essentially a two-point autocorrelation function 
normalised by the variance. 

The other statistic used is an rms difference statistic: 



a 2 = 



ST 
~~T 



ST 
~T 



-1 2v 



(10) 



using the same notation. 

Both statistics are means for all six circle pairs for the 
given point in parameter space (Z , b , 9, ±1, a). 

For genuine matched circles, S should be larger than 
for surrounding regions in parameter space, while a should 
be smaller than for surrounding regions. 

4. Results 

4.1. Search 

The full parameter space [eq. (JSJ] is searched with a res- 
olution of 2°. Since the parameter space is large, not all 
values of S and a are stored. Instead, over the 9 dimen- 
sion, values of S and a are stored for the two values of 
9 which maximise S and minimise a respectively. These 
two values of 9 may be, but are not necessarily, identical. 
Cornish et al. (2004) write S'max to indicate that S is a 
maximum value, not a mean value, but this is not strictly 
identical to the value used here, since values of both S 
and a are stored for both the cases of maximum S and 
minimum a. 

Taking the same approach as Cornish et al. (2004) , the 
values of S and er are plotted as a function of a, for the 
Clifford translations of both orientations, i.e. rotations of 
±36° and for a zero rotation. 

These values are shown in Figs0HHl 

There is clearly a cluster of points at ~ 10° which 
show an excessly high values of S (Fig. and excessively 
low values of a (Fig. [TJ for the left-handed rotation, in 
contrast to the general distributions of these statistics in 
the right-handed and zero (unphysical) cases. 

The general increase in the value of S as a decreases 
is similar to its behaviour in the figures of Cornish et al. 
(2004), and consistently, there is a corresponding decrease 
in a values as a decreases. 

Could this cluster of points indicate a real signal? 

If they correspond to a real signal, spread out over 
about a degree due to smoothing of the ILC map, to origin 
of the signal from the naive Sachs- Wolfe effect at super- 
degree scales, and possibly due to intrisic variations in Q\ 
of Afi A ~ 6xl0" 5 [eq. ©], then 

(i) they should represent adjacent points in parameter 
space (Z^, 9, ±1, a), and 

(ii) each of the six individual circle pairs should show ex- 
cess correlations and low differences near the optimal 
point. 
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angular radius of circle (degrees) 

Fig. 4. Values of the correlation statistic S [eq. @] as a 
function of matched circle angular radius a for each pair of six 
matched circles implied by an orientation of the fundamental 
dodecahedron, where the rotation of the Clifford translation 
(screw motion) mapping from a circle to its image is —36° (left- 
handed). The values of a are discrete, but a small random offset 
in a filling the 0.5° interval between discrete values is used for 
to better display individual points. The intervals in l^,b^ and 
9 are 2°. Points with excess correlation appear at ~ 10°. 



m o 




angular radius of circle (degrees) 

Fig. 5. As for Fig. 0] for a right-handed rotation. If a signal 
is present, it should be present for a left-handed rotation or a 
right-handed rotation, but not for both. 



4.2. (i) location of the ~ 10° excess 

The twenty points with the highest values of <S* in the range 
5° < cv < 15° all lie in one of the two ranges 

Z 11 = 253 ±4° 
6 11 = 64 ± 1° 
6 = 57.0° 

(the point with the lowest S value is at 59.0°) 



m o 




angular radius of circle (degrees) 

Fig. 6. As for Fig. 2] for a zero rotation. This is unphysical, 
so only noise signals should be present. 



or 



l U = 50±3.5 C 
b U = 51.0° 
9 = 25 ±2°. 



(12) 



(11) 



These two triples (I , b ,9) correspond to the same 
orientation of the dodecahedron — they are two face cen- 
tres for a single orientation. As mentioned above, the full 
parameter space is covered a little more than once, in or- 
der to be sure that all orientations are tested. 

The difference statistic a varies more rapidly with a, 
so a narrower range around ~ 10° is required to see if a 
isolates similar points to S. 

Of the twenty points with the lowest values of a in the 
range 8° < a < 15°, all but four lie in the ranges of l|llfl 
and (|12f) . except for four points with 

l U = 57.5° 
b 11 = 80.0 ± 1° 
6 = 67 ±2°. (13) 

However, the four points of (|13|l have correlations 0.27 < 
S < 0.30, much lower than for the other sixteen points, 
which have 0.50 < S < 0.61. This solution is not a face 
centre of the same dodecahedron as (|11|) and (|12fl . 

Clearly, only the solution Ijllj l and (|12H is simultane- 
ously indicated by both statistics. 

4.3. (ii) individual circle pairs 

If the ~ 10° signal is really due to the Poincare dodecahe- 
dral topology, then it should appear independently in each 
of the six circle pairs. If it is due to non-Gaussian and/or 
non-Poissonian properties of the ILC, it is unlikely to ap- 
pear in all six pairs. 

To test this, the correlation statistic S is calculated 
separately for each of the six individual pairs within the 
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angular radius of circle (degrees) 

Fig. 7. As for Fig. [I] showing values of the difference statistic 
a [eq. ill It ] as a function of matched circle angular radius a for 
each pair of six matched circles implied by an orientation of the 
fundamental dodecahedron, again for a left-handed rotation. 
A cluster of points with smaller differences a than the general 
distribution is visible at ~ 10°. 




10 20 30 40 

angular radius of circle (degrees) 
Fig. 8. As for Fig. for a right-handed rotation. 

Table 1. Sky positions of the six face centres for the do- 
decahedron which shows excess values of the correlation 
statistic S. A face number i, galactic longitude, latitude 
and estimated circle radius a (all in degrees) are listed. 
The other 6 faces are directly opposite. The orientation of 
the 36° screw motion between faces is left-handed. 




10 20 30 40 

angular radius of circle (degrees) 
Fig. 9. As for Fig. [7] for a zero rotation. 



range 

I 11 = 253 ±5° 
b 11 = 65 ±1° 
6 = 57.0 ±1°, 



(14) 
(15) 
(16) 



with a 0.5° resolution. 

Since the domain in most parameters is small, the 
mode for any given a value is plotted rather than a scatter 
plot of the individual values. 

Figs HOlfTI and IT51 clearly show that the peak correla- 
tion at around 10° is present in all six circle pairs for the 
left-handed screw motion. 

Correlations for a few individual peaks reach similar 
values for the right-handed and zero (unphysical) screw 
motions, but not for all six pairs simultaneously. 

The maximum of the individual curves in Fig. I1UI is 
used to estimate the circle radius. 

As shown in Tabled these radii all agree to within one 
degree, i.e. a = 11 ± 1°. 

4.4. Matched circles 

What are the actual values of the temperature fluctuations 
along the identified circles? The identified circles for the 
central values listed in eq. (|16|) are shown in Figs ll3hH%l 

These clearly show that temperature values around 
these circles are highly correlated. 



i 


,11 ■ o 

I in 


,11 • o 

o m 


o 

a m 


1 


252.4 


64.7 


9.8 


2 


50.6 


50.8 


10.7 


3 


143.8 


37.8 


10.7 


4 


207.5 


9.5 


10.7 


5 


271.0 


2.7 


11.8 


6 


332.8 


25.0 


10.7 



4.5. Point source catalogue 

If the matched circles are due to topology, then any fore- 
ground point sources close to them in celestial position 
should be present in one circle but not the other. 

Of the 208 point sources listed in the WMAP point 
source catalogue (Bennett et al. 2003), only one is 
within 0.5° of the matched circle solution. This source 
is GB6 J1635+3808 at (l U = 61.1°, 6 11 = 42.3°) and is 
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angular radius of circle (degrees) 

Fig. 10. Values of the correlation statistic S calculated for 
the six individual pairs of circles for dodecahedrons with ori- 
entations in the apparently favoured range [eq.lJTBJ], shown as 
the mode at each circle radius a. This is for a left-handed screw 
motion. All six pairs have high correlations at around ~ 10°. 



located within 0.1° of the circle shown by the solid curve 
in Fig. [H 

The point source appears to be present in the curve (at 
position along the circle of about lOOO/i -1 Mpc, shown by 
a pair of overlapping circles). A removal of flux in order to 
correct for the presence of this point source would improve 
the match between the circles. 

5. Discussion and Conclusions 

It is clearly premature to claim a highly significant de- 
tection of the topology of the Universe based on just one 
simple analysis of the first year WMAP ILC cosmic mi- 
crowave background map. 

However, the plots are striking and it seems prudent 
to release them to the scientific community while a com- 
panion paper is prepared with formal statistical analyses. 

Whether or not the matched circles found are just coin- 
cidence or due to global geometry, it is clear that tempera- 
ture fluctuations around 12 dodecahedrally spaced circles 
of radius 11 ± 1° in the WMAP ILC map correlate unusu- 
ally well in their respective pairs when a phase shift of 36° , 
corresponding to a left-handed screw motion, is applied. 

If simulations can show that this is a fairly likely occur- 
rence due to Gaussian fluctuations in an infinite flat uni- 
verse, then this would show that genuine matched circles 
will be even harder to distinguish from spurious detections 
than was previously thought. 

However, Vielva et al. (2003) have found non-Gaussian 
fluctuations at about 10° - extremely close to radius of 
the matched circle radius. While it is not obvious how 
this non-Gaussianity detection should relate to the radius 
of matched circles, since correlations are between circles, 
not along individual circles, this would also complicate 
simulations, since they would need to be consistent with 
observational analyses like this one. 




angular radius of circle (degrees) 

Fig. 11. As for Fig. EJD for a right-handed screw motion. 
Although correlations of some circle pairs are high, there is no 
case where all six are high. 




angular radius of circle (degrees) 

Fig. 12. As for Fig. 1101 for a zero (unphysical) screw motion. 
Although correlations of some circle pairs are high, there is no 
case where all six are high. 



Cornish et al. (2004) avoided circles of radii smaller 
than 25° because of the risk of false positives, and car- 
ried out extensive simulations of what signal would be 
expected from a genuinely multiply connected universe. 

However, these simulations risk the problem of cos- 
mic variance and logical circularity. If the Universe really 
is detectably multiply connected, then it would be fairly 
reasonable that the density perturbations (eigenmodes) 
on the largest scales are somewhat affected by the global 
physics of the Universe, so that modelling them on the 
basis of Gaussianity and random phases may lead to sta- 
tistical statements which are incorrect, because they talk 
about an ensemble of likely universes with different sta- 
tistical properties to the real Universe. 

The cosmic variance problem is that our actual 
Universe is just one realisation — the number of pertur- 
bations on large scales are too small for the large number 
theorem to make statistics of ensembles valid. 
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Fig. 13. Temperature fluctuations in the WMAP ILC data around an 11° radius circle centred at position i = 1 in Tabled 
(continuous curve) and its opposite circle (dot-dashed circle), matched by a screw motion with a left-handed rotation of 36° 
shown against the distance around each circle for (flm = 0.3, S7a = 0.7). These distances are a slight underestimate of the true 
distances for the hypothesis, since Euclidean calculation is used for a 'slightly' positively curved space. The two curves shown 
are 8T/T In the general analysis in this paper and in this and following plots, point pairs where either point lies at a galactic 
latitudes with |fe^| < 2° and or within 20° of the galactic centre are excluded. 



An apparently strong heuristic argument against these 
matches being physical is the need for extreme fine-tuning. 
If the claimed matched circles are due to topology, then 
the in-radius of the Universe happens by chance to be 
about cos(ll°) « 98% of the distance to the surface of 
last scattering. Intuitively, this is difficult to accept. 

However, as pointed out by Luminet et al. (2003), 
in the case of a positively curved Universe, especially 
if fltot ~ 1-01, there is necessarily a fine-tuning, since 
the curvature scale is only about three times the matter- 
dominated horizon size. 

Moreover, the non-zero cosmological constant, 
0.7, is well established observationally (though see 
Blanchard et al. (2003); Vauclair et al. (2003) for a minor- 
ity viewpoint) and definitely requires fine-tuning of some 
sort. 

Attempts have been made to link a non-zero cosmo- 
logical constant with detectable cosmic topology, but so 
far no obvious successes have been found. 

Independently of theoretical arguments, the results 
presented in this paper are testable by several observation- 
ally based methods which do not require assumptions on 
hypothetical statistical ensembles of universes, and each 
should potentially be able to improve the signal if it is 
cosmological in origin: 

(i) an attempt to separate out the naive Sachs- Wolfe ef- 
fect, the doppler component, and the integrate Sachs- 
Wolfe effect, and analysis at higher resolution 



(ii) improved removal of foregrounds, either indepen- 
dently of the hypothesis, or using the matched circles 
hypothesis to predict foregrounds which intervene in 
one circle but not the other (cf. Sect. 14. 5|) 

(iii) polarisation data from Planck. 

Similarly to the prediction by Luminet et al. (2003) 
that l Q tot ~ 1-013 > 1', Fig. □ shows that for Q m = 
0.28 ± 0.02, the total density parameter must be f^tot ~ 
1.010 ± 0.001 for these matched circles of radius 11 ± 1° 
to be cosmological in origin. A larger uncertainty in O m 
would correspondingly increase the uncertainty in fltot: 
but the prediction that iltot > 1 remains: the fundamental 
dodecahedron of the Poincare dodecahedral manifold is 
that of a manifold of positive curvature. 
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Fig. 14. As for Fig. 1131 for position i = 2 in Tabled This is the only one of the six circle pairs for which one of the Bennett 
et al. (2003) point sources lies within 0.5°. This source is GB6 J1635+3808 at (I ^ — 61.1°, = 42.3°) and is shown by a pair 
of overlapping circles at (arbitrarily) ST/T = 0. Its sky position is within 0.1° of the circle shown by the solid curve, suggesting 
that it is a contaminant present in the solid curve (Sect. 14.51 . 
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Fig. 15. As for Fig. 1131 for position i 
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